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TRANSIENT RESPONSE FROM THE LYAPUNOV STABILITY EQUATION 


By 

William G. Vogt 


SUMMARY 


For a certain class of linear and nonlinear asymptotically stable systems 
and a certain class of Lyapunov functions v(x) , estimates of transient response, 

-v(x) v -v(x) 

= max [ — — 1 , y = min [ — — ] 

1 v(x) J v(x) 

X x 


A 

r 


are related to the real parts of the characteristic values of the system matrix, 
A, of the linear approximation system. It is shown that 


y 5 -2 min [Re X. (A) ] 
ji 

y S -2 max [Re A. (A) ] 

A V 

and a formula is given to choose v(x) such that y, y approach these limiting 
values as closely as desired. 


INTRODUCTION 


In stability studies of control systems describable by a vector 
differential equation of the type 

x = f(x) JL(0.) ~ 0. (1) 

where x and f^(x) are n -vectors and f(x) satisfies certain conditions [1] 
assuring the existence, uniqueness, and continuity of solutions to equation 
(1) , a common procedure is to obtain a Lyapunov function v(x) which, along 



with its derivative v(x) , satisfies conditions sufficient to assure stability or in- 
stability of the null solution x = 0^ of equation ( 1) . More specifically, for most 
control systems, asymptotic stability is desired. A sufficient condition for 
asymptotic stability of the null solution (NS) of equation (1) is the existence of 
a Lyapunov function v(x) which, along with its derivative v(x) , satisfies the 
following conditions in some neighborhood of the origin [2-4] : 


(a) v(0) = 0 

(b) v(x) > 0, x^|0 

(c) grad v(x) is continuous 

(d) v(x) = fgrad v(x) ]* i(x) < 0, x *0. 


If such a v(x) exists, the asymptotic stability of the NS of equation (1) 
is assured. For such v(x) , it has been shown that an estimate of the transient 
response of equation (1) can be obtained [2], This estimate consists of de- 

V A 

termining two constants y, y by the following formulas 


v 

y = mm 
x 


-v(x) 

v (x) I; 


A 

y = max 
x 


-v(x) 

v(x) 


(3) 


in some region 0< | x | <r,r>0 where the conditions of equation (2) are 
satisfied. Denoting by _0( t; x 0 ) a solution of equation ( 1) starting from x 0 
at time t=0, the time behavior of v(£(t; x 0 ))must satisfy [2] 

v(x 0 ) exp [ - yt] g v(0(t;x o )) g v(x 0 ) exp [-yt]. (4) 

A V 

Thus y and y correspond in some sense to the largest and smallest time 
constants of the system described by equation (1) , but the exact relationship 
is not clearly understood. 

A V 

The significance of y and y is obscured largely because the Lyapunov 
functions are not unique. Different Lyapunov functions yield different values 

for y and y. These estimates depend on the system of equations and the 
Lyapunov function chosen. 


The object of this paper is to relate the estimates in question to a 
certain class of equation (i) and a certain class of Lyapunov functions 
satisfying equation (2) . 
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FORMULATION 


System 

The following investigation is limited to that class of systems of 
equation (1) for which the right sidej[(x) satisfies the following conditions 

f(x) = A x +jg (x) (5) 

Re[X.(A)]<0 1=1,2, . . .n. (6) 

In words, _f(x) is limited to that class of vector functions, each com- 
ponent of which is expressible in a convergent power series expansion in some 
neighborhood of the origin and the linear approximation of which has only neg- 
ative real part characteristic values. If A is a matrix that satisfies equation 
(6) it is called a stability matrix. 

The ordinary linear constant coefficient differential equation which 
approximates equation (5) in some sufficiently small neighborhood of the NS 
of equation (5) is given by 

x = Ax (7) 

and the solutions are designated by 0 (t; xq) • 

Lyapunov Functions 

Attention is restricted to Lyapunov functions of the form 

v(x) = x'Vx +vj(x) (8) 


where V is a real, symmetric, positive -definite (RSPD) matrix and Vi(x) 
contains only terms of degree greater than two in the variables x t , X 2 , . . . .x^. 

Thus , the first approximation to v (x) is given by 

V(x) = x' V x . (9) 

The derivative of v(x) along solutions to equation (1) is given by 



where vj(x) consists only of terms of degree greater than two in the variables 
x 1# x 2 , . . . , x^. The first approximation of v(x) is therefore given by 

V(x)=x'(A'V + VA)x . (li) 

For v(x) , v(x) to satisfy the requirements of equation (2) in some 
neighborhood of the origin, V(x) must be negative definite, thus implying 
that 


A' V + V A = -2R (12) 

where R is an RSPD matrix. 

Equation (12) is referred to as the Lyapunov stability equation. To 
illustrate its fundamental importance, a theorem is stated. 

Theorem 1 [2] : A is a stability matrix if, and only if, V, as a solution to 
equation (12) , is unique and RSPD wherever R is RSPD. 

Thus, instead of choosing an RSPD V in equation (9) such that the 
resulting R in equation (12) is RSPD, choosing any RSPD R and solving for 
V in equation (12) settles the stability question of the NS of equation (7). 

If there is asymptotic stability of the NS of equation (7) , V provides an 
appropriate Lyapunov function V(x) in equation (9) which meets all the con- 
ditions of equation (2) . Thus, in the system formulation, equation (6) can 
be replaced by the condition that V, as a solution to equation ( 12), is unique 
and RSPD wherever R is RSPD. 

Note that the conditions of equation (5) and (6) are sufficient to 
assure the existence of Lyapunov functions of the type required by equation 
(8) (3]. With this formulation, the problem originally posed can be reduced 
to a simpler problem. 


REDUCTION OF THE PROBLEM 


A V 

The original problem was to determine how the constants y and y 
defined in equation (3) are related to the systems and the Lyapunov functions 
chosen. From the definition 


y(x) 


-v(x) 

v(x) 


(13) 
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( 14 ) 


it is clear that for 0 < j x J < r 

A 

y = max [y (x) ] 
x 


v 

y 


min [y (x) ] . 

x 


Substitution of equations (8) and 


V(x) + v t (x) 

y(x) = _ 

V(x) + Vi(x) 


( 10) into equation (13) yields 
v t (x) 


1 + 


V(x) 


= 2T(x) 


1 + 


vi (x) 
V (x) 


where 


r(x) = 


1 V(x) 

2 V(x) 


x ' R x 
x'Vx 


(15) 


(16) 


From equation ( 15) and the nature of the Lyapunov function and its derivative 
in equations (8) and (10) , it is clear that for all 0 < | x | < r* < r there 
exists an e (r^) such that 


2T (x) (l-2e) S y (x) S 2F (x) (l+2e) (17) 

and further that 



lim e (ri) 

= 0 as rj 0. 

Thus, 

for 0 < | x 

|<r 


A A 

V v 


y = 2 r ; 

y S 2T 

where 

A 

x ' R x v 


T = max 

[ x'Vx 1 ; r - mm [ 


X 

— X 


(18) 


(19) 


( 20 ) 


Clearly, for the class of systems and the class of Lyapunov functions 

considered, r and ^ give the limiting values of $ and y. The original problem 
is thereby reduced to the consideration of the following: 
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Problem: Given a stability matrix, A, what are the interrelationships between 

a v 

A, R, and V satisfying equation (12 ) in terms of r and T defined by equation 
( 20 )? 


The answer to this reduced problem, while still incomplete, is 
scattered throughout many publications [5-7] , often in forms which are not 
suitable to the problem described here. The treatment below is an attempt 
to organize this material to give an answer to this problem and at the same 
time to introduce some new results. 

a v 

Note that equation (4) is not satisfied with 2T and 2r substituted for 

A V 

y and y. The relations of equation (19) indicate that the fastest (slowest) 
system "transient" in terms of equation (4) is as fast (slow) or faster 

A V 

(slower) than that characterized by exp [-2r t] (exp [ -2F t ]) and further 
that near the origin the fastest and slowest response are characterized by 

A .V 

"time constants" l/2r and 1/217 respectively. 

PENCILS OF QUADRATIC FORMS [5] 

Definitions and Theorems 


The real form 

x'Rx-x’Vx (R = R', V = V r ) (21) 

is called a pencil of quadratic forms. If x' V x is positive definite, the 
pencil is called regular. The equation 

|R - XV | = 0 (22) 

is called the characteristic equation of the pencil. If A^ satisfies equation (22) , 

it is called a characteristic value of the pencil. If there exists a real vector, 
z., such that 

R z. = A. Vz. , (23) 

- -i l 1 

z. is called a characteristic vector belonging to the pencil. The following 
- 1 

theorem holds. 
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Theorem 2: The characteristic equation (22) of the regular pencil of forms, 
equation (2i) , always has n real roots, A., with corresponding real character- 
istic vectors satisfying equation (23) . The z can be chosen such that 

z\ Yz = 5., (i,k = i, 2 n). (24) 

1 1C IK 

Since [ v| # 0, the matrix R V has the same characteristic values as 
the pencil, equation (2i), does. Hence, these characteristic values will be 

designated by A. (R V -1 ) , i = 1, 2, . . . , n. The matrix R V -1 has n real roots 

and is similar to a diagonal matrix [5] . 

Let Z designate the matrix composed of the columns^, which satisfy 
equation (24) . Then 


Z' V Z = E (Identity Matrix) 

(25) 

Z'RZ = diag { A.(R V -1 )} 

(26) 


Extremal Properties 

Order the characteristic values of R V -1 as 

min [A. (R V -1 )) = A^RV -1 ) S... ^ A ( R V' 1 ) = max [A. (R V" 1 ) ] 
i i 


and let the corresponding characteristic vectors satisfying equation (24) be 
Zj, _z 2 z . With this ordering, the following holds [5], 


x' R x 

A 1 (RV~ 1 )= min [ x , v “ “1 
x — 


z[ R z t 
zi Y z t 


A (R V -1 ) 
n 


max 

x 


x ! R x 
r x' Vx 


z' R z 
-n n 


z' V z 
— n n 


(27) 


(28) 
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A V 

Relations to f and f 


From the definitions of T and r given by equation (20) it is clear that 


A 

r 


X (RV _1 ) 
n 


(29) 


v 

r 


h (5 Y' 1 )- 


(30) 


Thus the problem has been reduced to finding the minimum and maximum 
characteristic values of the matrix R V -1 where R and V are related to A by 
equation (12). 

Remark : In this development, it was not assumed that R was an RSPD matrix, 
but merely that V was RSPD. Thus if R is chosen RSPSD (real, symmetric, 
positive -semidefinite) such that V in equation (12) is RSPD, the same relations 
hold. 


RELATIONS BETWEEN V (RV)" 1 AND X.(a) 

i — J 


Notation 


Adopt the notation 

A.(A) = ffl.(A) +1 /3. (A) (i_ = aTTI). (31) 

Note that a, < 0 since A is a stability matrix. Order the eigenvalues of A as 
follows: 


max [a. (A) ] = ct^ (A) ^ g a. (A) = min [of. (A) ] . (32) 

i i 

For all A. (A) having a. as a real part, some of these may have simple 

elementary divisors and some may not. In the case that there is at least one 
non -simple elementary divisor associated with the roots having oq( A) (a^(A) ) 

as a real part, designate one of these as ^(A) (^ fl (A) ) . Then «i(A) (« n (A)) 

either corresponds to characteristic values having simple or non -simple ele- 
mentary divisors depending on whether A has or has not at least one non -simple 
elementary divisor corresponding to roots with real part 04(0^) . 
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Complex-Valued Vectors 

For convenience in the later development, the following theorem is 


stated. 


Theorem 3: If H is any real symmetric matrix and w is any complex -valued 
vector, a real vector x can always be found such that 


x'Hx = w* Hw (w* = w') . (33) 

Proof: Let 

D(x, w) = x'Hx - w* Hw (34) 


and w =u+iv, x = u + b v where b is a real constant which is to be selected 
such that D(x,w) = 0. Note that w* Hw is a hermitian form and thereby takes 
on only real values. 

D(x,w) = (u + bv)'H(u + bv) - (u -_i v) ' H(u + _i v) 

= u ' H u + 2b v ' H u + b 2 v ' Hv-u'Hu-v'Hv 
= b 2 v ' H v + 2 v ' H u b -v'Hv 


H b is selected as follows, D(x,w) will be zero, 
. 0 if v’ H v = 0 

-v ' H u + [ (v ’ H u) 2 + (v ' H v) 2 ] 


' if v* H v #0 


(35) 


Choosing b according to equation ( 35) makes D(x,w) = 0. Thus x = u + b v is 
the real vector such that equation (33) holds. 


Relations Between x.(BY) 1 And x.(A) 

i j 

Theorem 4: Given a real stability matrix A, and any RSPD or RSPSD R and 
RSPD V satisfying equation ( 12) , the following inequalities hold: 

MR V" 1 ) ^ - «i(A) = -Re [^(A)] (36) 
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( 37 ) 


A (R V" 1 ) 5 -oi (A) = - Re [A (A) ] . 
n n — n — 1 


Proof: By equations (27) and (28), and Theorem 3, it is necessary only to 
show the existence of complex -valued vectors w. such that 


- a 

l 


w * Rw. 

— l 1 

w* V w. 

— l 1 


(i = 1, 


n). 


(38) 


From equation (12) 


w* (A'V + VA) w. = - 2 w.* Rw. . (39) 

— l — — — l ~ l — l 

Let w. be a characteristic vector of A; i.e. A w. = A. (A) w.. Then equation 

- l — ii— — 1 

(39) yields, since A* = A* , 


(Aw.)* Vw. +w.* V (Aw.) = -2 w* Rw. 

1 7 1 -i i -i 1 


(A. + A.) w.* V w. = -2 w* R w. . 

l i—i 1 —i 1 


(40) 


Since A. + A. = 2 a.. Theorem 3 shows that there exists a real vector x. 

... i 1 i _1 

with 


- Oi. 
i 


x : r x 


x! Vx. 

— l 1 


i = 1, 2, 


n. 


(41) 


Since equation (41) holds in particular for i = 1, n and since the x. may not be 
the principal vectors z^, equations (36) and (37) must be true. Thus the theorem 
is true. 


Can the equalities of equations (36) and (37) be made strict for some 
choices of R and V? The answer to this question must await the results of 
the next section. 
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A TRANSFORMATION OF THE LYAPUNOV STABILITY EQUATION 


If a conjunctive transformation by the complex matrix T, | T | + 0 
is performed on equation (12 ) yielding 

T* (A'V+VA) T = -2 T* R T, (42) 

equation (42) may be rewritten as 

(T _1 A T)* T* V T + T* V T (T -1 AT) = -2 T* R T 


or 


B* U + UB = -2Q. (43) 

Where 

B = T _1 A T, U = T* V T , Q = p RT, 

it is clear that the definiteness properties of U and § are the same as those 
of V and R, respectively, and the elementary divisors of B are the same as 
those of A. 


In particular, let T = T(e) transform A to its Jordan normal form 


T" 1 A T = diag {J 4 J } = J (e) 

s 


J. 

1 


X.(A) e 0 

. ’ € 

0 X. (A) 


X. (A) E. + eH. 

l ~ —i — l 


(44) 


(45) 


where E. and H. are of dimension p. x p. corresponding to the degree of the 
elementary divisor of A,(X. - X.(A))^i. 

In equation (43) , choose Q to be Q = diag (Qi, Q } 

s 

where 


S.(€) = - Q'.(A) E.+ (|)(W + H i )- 
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Thus, U = E (the identity matrix). 


For e > 0, but sufficiently small, it is clear that Q so defined is 
positive -definite. The characteristic values of the pencil (Q - AE) are 
the roots of the equation 


j Q - AE | = 0 


which are the roots of the equivalent expression 
s 


n 

i=i 


AE. - Q. = 0. 

—i —i 


( 46 ) 


(47) 


By a theorem of Gershgorin [8] , the roots of each of these factors of 
equation (47) must be within an e distance of the corresponding - A. (A) ; i.e. 


-o'. (A) -e ^ A. (Q.) S - “.(A) + e (j. = 1, 2, . . .p.) . (48) 

In particular, for a J. corresponding to a simple elementary divisor, 

A. (Q.) = -a i(A). (49) 

Ji 1 

Since R V -1 is similar to Q, i.e. 

R V" 1 = T* _1 Q T _1 T = T* _I Q T* (50) 

the A.(R V -1 ) are the same as the A.(Q) . The foregoing, along with Theorem 
J 1 

4, has demonstrated the following theorem. 

Theorem 5: Given a stability matrix A, and an e > 0, it is always possible to 
find an R and V satisfying equation (12) such that 


Q-i - e g At(R V -1 ) ^ - a t 

(51) 

a = A (R V" 1 ) ^ - a + e. 

(52) 


n n n 


In particular, choosing 



fp — 1 'pi— l ip-lj 


(53) 
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where T transforms A to J(e) as in equation (44) , yields the required R and V. 
Note that R and V depend on e . Consideration of equation (49) yields the 
following corollary. 

Corollary 1 : Given a stability matrix A, an RSPD R and V satisfying 
equation ( 12 ) can be found such that 


A i (RV- 1 )= -a t (A) 

(54) 

[ V5Y' 1 ) =-« n ( A)] 

(55) 


(a) if and (b) only if the elementary divisors of A corresponding to the A j(A) having real 
part aq(A) [a (A) ] are all simple. 

Proof: If they are all simple, the implication (equation (54)) follows directly 
from equation (49) , establishing (a) of the corollary. To establish (b), it is 
only necessary to show that if one A. (A) with Re [A^(A) ] = o; 1 (A) has a non- 
simple elementary divisor, there is at least one vector, say x, such that 

x'Rx 

jryr < . 




Let u and v be complex such that 

Au = [ 04 (A) + _i j3.(A)] u - A.(A) u (56) 

A v = [ct^ (A) +_ii3.(A) ] v + u = A.(A)v+u. (57) 

Thus A. (A) has a non-simple elementary divisor and real part q' 1 (A). Let 
w = u + e v where e is real. From equation (12) 


w* (A' V + V A) w = -2 w* R w . 

From equations (56), (57), and (58) 

[ Ai + A.] w* V w + 2 ew* V w - 2 e 2 (Re[v* V w] ) = -2 w* R w 
or 

Re [v* V w ] w* R w 

- a t - e (1 - e ^ y w ) = w* V w * 


(58) 


( 59 ) 
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For sufficiently small e, say |e<)| > | e | > 0, the parenthesized term on the left 
of equation (59) will be positive. Thus, for |eo| > € > 0> 

w* R w 

w* Vw < ~ “1 • 


By Theorem 3 , a corresponding inequality will hold for some real x. This 
implies 


Ai(R V" 1 ) < - a>i 


(60) 


which completes the proof of Corollary 1. 

OTHER BOUNDS ON THE X.(B¥' 1 ) 

I 

Bounds on x (RV 1 ) 

1 

The preceding sections have shown that Aj(R V -1 ) is bounded above by 
- Qij( A) . If A is a stability matrix and R and V are to be RSPD matrices re- 
lated by equation (12) , Aj(R V -1 ) must be bounded below by zero. Thus 


0 < Ai(R V -1 ) < - 04(A) 


( 61 ) 


If the conditions are relaxed to include RSPSD R, while maintaining the re- 
quirement that V be unique and RSPD as a solution to equation ( 12) , equation 
(61) takes the form 

0 ^ MRV- 1 ) g-Q! 1 (A) . (62) 

According to a theorem by LaSalle f9] , asymptotic stability of the null so- 
lution of equation (1) is assured if the conditions of equation (2) hold with 
(d) relaxed to 


(di) v(x) SO 


(63) 
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provided the set of x such that v(x) = 0 contains as a largest invariant set of 
equation (i) , the null solution, x = 0. In terms of V(x) and V(x) , this re- 
quires that the null space of R contain no invariant space of A except x = 0. 
Provided these conditions are met, R can be chosen RSPSD, justifying the 
equality on the left of equation (62) . 

Bounds on X (RV _I ) 
n 

A (R V -1 ) is bounded below by - a (A) . It is also bounded above as 
n n — 

long as R is at least RSPSD. In particular , 

n 

-« n (A) 5A n (R V" 1 ) <-tr (A) = - £ o'. (A) (64) 

i=i 

Equation ( 64) follows from equation ( 12) by considering 

A’ + V A V -1 = - 2(R V -1 ) 

tr (A' + VA V _I ) = 2 tr (A) = - 2 tr (R V -1 ) 
or 

- tr (A) = tr (R V" 1 ) . 

n 

Since X. (R V -1 ) > 0 and since tr (R V -1 ) = ^ X. (R V -1 ) , equation (64) 

1 i=i 1 

follows directly. 

RELATIONS TO T AND f 


The following estimates have been obtained. 

-« n (A) i r g-tr(A) (65) 

o S r < - o-^A). (66) 
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These estimates are valid for any RSPD R and V which satisfy equation ( 12) . 

a v 

If a judicious choice of R and V are made, T and T can be made as 
close to - Q!^(A) and - 04(A) as desired, i.e. for e > 0 . 

A 

- (A) % r S - a n (A) + e (67) 

- 04(A) - e ^ T $ - 04 (A) (68) 

Such a V is given by equation (53) . 

In equation (67) , e can be set to zero if all A. (A) with Re [A.(A) ] = 

a (A) have simple elementary divisors; similarly, for e in equation (68). 
n “ 


RELATIONS TO r AND r 

From equation (19) it is clear that for the class of systems and the 
class of Lyapunov functions, 

A A 


y > 2 T > - 2 o; (A) 

(69) 

v v 

0 Sy <2T S - 2 04(A) . 

(70) 


Equations (69) and (70) can be thought of in two ways. If y and y are 

known, they are estimates of - 2 a (A) and - 2 04(A). Conversely, if -a (A) 

n a n 

and -04(A) are known, they serve as limiting values for y/2 and for 7/2 for all 
Lyapunov functions of the class considered. 


EXAMPLES 


To illustrate the theory just stated, several examples are given. 
For ease in calculation and clarity in presentation, only the most simple 
matrices are selected for A. 
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The Lyapunov stability equation is 


A' V + V A = - 2 R. 


(12) 


In each example below, for the A and R given, V is solved for and the A^V 1 ) 
are calculated. 

Example 1: 


A = 


-a i 

0 


0 

-a 2 




1/ai 


V 


2k/ (aj + a 2 ) 


2k/ (at + a 2 ) 
r 2 /a 2 


where a 2 ^ a^ > 0, r 2 - k 2 > 0. After some manipulations 


Xi(Rr‘) = a i -K(a 2 -ai) ^ a t 
A 2 (R V *) = a 2 + K(a 2 -ai) = a 2 


where 


K 


i [ 


r 2 -k 2 


4aja 2 

(at+a 2 ) 2 


1 

2 


- 1 ]. 


Observe that K = 0, if k = 0; K > 0 if r 2 - k 2 > 0, and K(a 2 -aj) = a^ if 
r 2 -k 2 = 0. 


Note that even if r 2 - k 2 = 0, that is, R is RSPSD, V is still RSPD and 
according to the theorem by LaSalle, this is sufficient for the asymptotic 
stability of x = 0^ of equation (7) , because the null space of R contains no 
invariant subspace of A. Also note that if r 2 - k 2 = 0, 
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AiCRV- 1 ) = 0 


Ag(R V ^) = + &2 


or, in other words, the limiting values on the A. (R V' 1 ) predicted by 
equations (62) and (64) are actually attained for this choice of R. 


Example 2: 



V = 


i/a 


(2ak+l)/2a 2 


L (2ak+i)/2a 2 (2a 2 k z +2ak+i)/2a 3 


where a > 0, r 2 -k 2 > 0. The A.(R V *) are 

1 

Ai(R V *) = a [i - (4a 2 (r 2 -k 2 ) + 1) 2 ] < a 

i - ~ 

A 2 (R v" ) = a [1 + (4a 2 (r 2 -k 2 ) +1) 3 ] > a 

It is clear that no value of k will make Aj(R V -1 ) = a = A 2 (R V -1 ) . 
However, given any e > 0, the relation 

a > A t (R V -1 ) > a - e 

a < A 2 (R V- 1 ) < a + e 
can be made to hold by choosing 
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] 


r 2 - k 2 > max [ 0 , — 

4a 2 e 2 

which illustrates Theorem 5 and its Corollary 1 . 

Also note that if r 2 -k 2 = 0, that is, R is RSPD, V is still RSPD, and since the 
null space of R contains no invariant subspace of A, this is still sufficient to 
show that x = 0 of equation (7) is asymptotically stable. In this case 


^(Rr 1 ) = o 

X 2 (R V- 1 ) = 2a 


(r 2 -k 2 = 0) 


which also illustrates equations (62) and (64). 
Example 3: 


A = 


-a i 

a 2 “ 


“1 

k “ 



R = 



^ a 2 

_a l _ 


k 

1 


V = 


a l + a 2 


2 2 
ai - a A a 2 k + a 2 


ajk 


a l 


a^k 

a^ + a^a 2 k + a 2 
a l 


where aj > 0, i-k 2 >0. The A.(R V -1 ) are 


^i(RV-‘) = a t [ 1 - ( 


2 i 2 

a 2 k 


a 2 + a 2 - a 2 k 2 


) 1 ^ a i 
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a 2 k 2 i 

X 2 (R V" 1 ) = at [ 1 + ( — ) ] > a t . 

a 2 + af - a 2 k 2 

Note that 

^(RV- 1 ) = aj = A 2 (RV _1 ) if k = 0 

and 

X,(RV J )= 0, \j(RV J ) = 2a, if k= 1. 


CONCLUSIONS 

A V 

It is evident that the estimates of transient response y and y will 
always give upper and lower bounds for the smallest and largest "time 
constants" of the system considered, but that such estimates may not be 
very good. Moreover, a judicious choice of the Lyapunov functions will 
give as good an estimate as desired. However , this choice of the required 
Lyapunov functions is intrinsically dependent on a knowledge of the charac - 
teristic vectors of the matrix A. This is clearly undesirable. Work is now 
underway to find some method to utilize digital computation to obtain this 
judicious choice of Lyapunov function without prior knowledge of the charac- 
teristic vectors of A. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Huntsville, Alabama, February 23, 1966. 
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